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Abstract* We suggest that cosmic rays, up to the highest energies obser- 
ved, originate in the galaxy and are accelerated in astrophysical shock 
waves. If there is a galactic wind, in analogy with the solar wind, we 
expect a hierarchy of shocks ranging from supernova shocks to the galac- 
tic wind termination shock. This leads to a consistent model in which 
most cosmic rays, up to perhaps 10*^ eV energy, are accelerated by super- 
nova shocks, but that particles with energies of 10 1 * eV and higher are 
accelerated at the termination shock of the galactic wind. 


Introduction. It appears that supernova blast waves are the major site 
for acceleration in the galaxy. While this mechanism is attractive, the 
energy to which particles can be accelerated is only of the order of 10 
eV (Lagage and Cesar sky, 1983). However, the observed cosmic-ray spec- 
trum extends well beyond this cutoff, and one is forced to consider 
separately the origin of very high energy particles. 

Two main factors limit the energy to which particles may be acce- 
lerated - their gyroradius (which is related to the mean-free path) and 
the time available. The time scale for particle acceleration is given by 
(e.g., Krymsky et al, 1979, Forman and Morfill, 1979, Axford, 1980) 
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where is the diffusion coefficient (which is a function of energy) and 
V fi h is the shock velocity. Clearly, since increases with energy, at 
any time there will be an energy above which the energy spectrum cuts 
off. We set the mean free path a constant of order unity times the gyro- 
radius in the ambient magnetic field. For particles of rigidity R (in 
electron volts) in a magnetic field of B (gauss) we have 

tacc ~ °[OQ 0,^ 

The velocity of a supernova shock wave may be written 


( 2 ) 


v sh 2 " < E sn/ M G> (3) 

where 1.33 for a pure blast wave and K 2.5 for evaporative solutions 
(Mckee and Ostriker, 1977), where E is the energy of the explosion and 
Mq is the mass of gas swept up by the blast. 

The value of Mq m ay be .found from the Sedov approximation. For a 
supernova of age t 5 times 1 (P years, and 1 , the maximum cosmic-ray 
energy in volts is J 

R(v) - 6.5x1013 83 <*> 

where fl 0 ~ Tl^XiaSn*, - E yj X 10*1*}$, and B - 3 *10 Bj Gauss. See also 
Lagage and Cesarsky(1983). The possible interact^n of cosmic rays with 
more than one shock does not appreciably alter this limit (Blandford and 
Ostriker, 1980). The spectrum of cosmic rays accelerated by supernovae, 
therefore is expected to steepen above a few times 10-*-^ eV. 

The observed cosmic-ray spectrum displayed in figure 1 shows that, 
although there is an apparent change in the spectrum in the vicinity of 
10^ eV, the acceleration process is still quite vigorous above this 
energy. It remains to find the source of the particles above approxi- 
mately 10*^ eV. 
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Figure 1. Spectrum 
of cosmic rays pre- 
sented by Linsley 
(1980)* 



The above discussion indicates roughly that shocks with a greater age 
and/or greater spatial extent can accelerate particles to a higher 
energy. Single supernovae shock waves presumably produce most cosmic 
rays, up to perhaps lO*^ eV. Super shells, containing the energy of about 
100 supernovae, all occurring in a short time, may also contribute. From 
equation (4) we find that the maximum cosmic-ray energy obtainable in 
supershells may be expected to increase only by a factor ~ 5 over that 
of a single supernova. 

The largest shock structure associated with our galaxy, and 
presumably the one with the greatest age, is the termination shock of the 
galactic wind. Presumably there will also be intermediate-scale shocks 
resulting from various processes in the galactic wind, similar in nature 
to those found in the solar wind. 

Acceleration at the Galactic Wind Termination Shock* The existence of a 
galactic wind flowing outward from our galaxy cannot be regarded as 
established, although it seems likely on the basis of several considera- 
tions* A mixture of gas, magnetic field and cosmic rays similar to the 
hot phase of the interstellar medium is expected to form a supersonic 
wind (e.g., Johnson and Axford, 1971, Ipavich, 1975, Jokipii,1976, 
Kopriva and Jokipii, 1983). In addition, the energy dependence of 
cosmic-ray composition suggests a general outflow from the galactic disk 
(see, «:.g. Jones, 1978, Jokipii and Higdon, 1979). We will assume the 
existence of a galactic wind in which the galactic-halo gas (the inter- 
cloud medium) flows outward, and terminates in a shock which can 
accellerate particles. We can estimate the properties of the galactic 
wind from some simple considerations. Irrespective of how the wind is 
driven, the energy source is the supernova explosions. We take a super- 
nova energy input rate into the interstellar medium of 

EgN = 10 ^ 2 ergs/sec 

and assume an energy redistribution between cosmic rays and galactic wind 
in the ratio 1:1 The energy input rate into cosmic rays, and hence into 
the galactic wind, is then of order 

E c * E w ■ MV w 2 /2 - 5xl0 40 ergs/ sec. 

We see immediately that the available time for acceleration at the 
galactic terminal shock is the lifetime of the galactic wind, which in 
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turn is T gal , the age of our galaxy, or 1.5 x 10 10 years (Hoyle 1961). 

I P - e . rgies of. gat tides accelerated in a. galactic wind terminal shock We 
US w$8 r P re V°us estimate for the energy input rate into the wind, E w = 5 
x IU ergs/sec. We express the mass loss rate M in units of the present 
diff us e galactic matter (5xlO^M Q ) divided by the lifetime, t ga i of the 
galaxy (1.5x10 yrs*). i.e., M 2xl0^^xF g/sec. where F is a scaling 
parameter. This yields a wind velocity V w ^5oo/Jf km/sec. A constraint is 
that V w > v E * 300 km/sec, the escape velocity. This implies that F < 4. 
in any case, mass loss should not be great, restricting F to perhaps < 1. 
Thus we do not expect mass balance problems. 

Because of galactic rotation, the field lines will be wound up much 
as m the solar wind. Hence, we use the Parker Archimedean spiral form 
ror the magnetic field and write for the magnetic field B at the shock in 
terms of that at the reference radius, B • 

® g Aj&'j o) 

wher e sec is the galactic rotation rate. Defining ^ and 

substituting for V w and setting t acQ = T g a ^ in equation (2) yields ) 

iv 7'^ (6) 

L i Tfjpn J 

In Table I are displayed R max for various using = 10 kpc, B n = 
^ ~ ** c * lo * c . e R 0 " 10 kpc is in agreement with a lower limit 

of 4 kpc for the static halo, using I.U.E. observations (Savage and 
DeBoer, 1981, Pettini and West, 1982). 

Table 1. Maximum rigidity (and energy of Fe nuclei) of cosmic 

rays accelerated at the galactic wind termination shock. 



Rghdcpc) 


E Fe (eV) 

0.2 

50 

7x10}£ 

1.8xl0 20 

0.1 

100 

3.2x10*8 

. 8 x 1020 

0.05 

200 

1.6x10*8 

.4x10 20 

0.02 

500 

6 .4x10*7 

.2x10 20 


The distance to the shock, R , should be about half the mean 
separation of the galaxies in our local cluster, or about 200 kpc. 

Finally, we note that particles accelerated at the terminus of the 
galactic wind will be subject to modulation by the galactic wind. Recent 
models of the solar problem (see e.g. Jokipii and Kopriva, 1979, Jokipii 
and Kota, 1983) suggest that the cosmic-ray intensity behaves roughly as 
if the particles were decelerated in a potential field equal to the 
electro-static potential difference between the heliospheric pole and 
equator,. Using our standard values and setting yields a value 
of 10 eV. This suggests that cosmic-ray particles with energies 
greater than this value will not be significantly modulated by the 
galactic wind, whereas lower-energy particles will be modulated. In this 
picture, galactic cosmic rays with energies less than about 10*4 eV are 
primarily- .accelerated by snpernova shocks, whereas those with energies 
above 10* J to perhaps 10 1 / eV are accelerated at the galactic wind 
termination. Those with intermediate energies may be accelerated by 
intermediate-scale shocks in the wind, in analogy to blast waves and co- 
rotating shocks observed in the solar wind. 

Further Consi derations and summary. The galactic disk cosmic rays repre- 
sent an average over many super-novae, at different stages (e.g., Bogdan 
and Volk 1983), whereas the galactic terminal shock is of more constant 
strength (and probably a higher Mach number), leading to a somewhat 
different source spectrum. 
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Since the acceleration process which we propose is quite slow, we 
must compare it with loss processes. For high-energy nuclei, interaction 
with the 2.7° K microwave background is the major consideration (Greisen, 
1966). Starlight is less important because of the distance from the 
galaxy. Puget, et al. (1976) showed that photodisintegration is the 
dominant process. Examination of their figs. 2 and 8 indicates ±hat 
protons and Fe have a loss time of t ea i for energy > 2 X 10 7 10 e y* 

Therefore, the present mechanism appears capable of accelerating cosmic 
rays to the highest observed energies. Finally, we suggest that 

the observed anisotropy in high-energy cosmic rays from north galactic 
latitudes may be related to the apparent motion of the galaxy, and the 
pushing of the termination shock closer to the galaxy in that direction. 

Hence, a single acceleration process in our galaxy may be 
responsible for the production of cosmic rays with energies up to 10* u 
eV. The generally uniform spectral characteristics of the cosmic rays 
over more than ten decades in energy are then easily understood to be a 
consequence of the same acceleration process - diffusive shock 
acceleration - and the fact that the shocks all share the same basic 
energy source - the supernovae explosions in the galactic disk. 
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